Josephson-like spin current in junctions composed of antiferromagnets and ferromagnets 
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We study Josephson-like junctions formed by materials with antiferromagnetic (AF) order parameters. As 
an antiferromagnet, we consider a two-band material in which a spin density wave (SDW) arises. This could 
be Fe-based pnictides in the temperature interval Tc < T < Tn, where Tc and Tn are the critical temperatures 
for the superconducting and antiferromagnetic transitions, respectively. The spin current jgp in AF/F/AF 
junctions with a ballistic ferromagnetic layer and in tunnel AF/I/AF junctions is calculated. It depends on the 
angle between the magnetization vectors in the AF leads in the same way as the Josephson current depends 
on the phase difference of the superconducting order parameters in S/I/S tunnel junctions. It turns out 
that in AF/F/AF junctions, two components of the SDW order parameter are induced in the F-layer. One 
of them oscillates in space with a short period f pj, ~ hvlJff while the other decays monotonously from the 
interfaces over a long distance of the order i, = Hvl2nT (where v, ^ and T axe the Fermi velocity, the 
exchange energy and the temperature, respectively; the subindex b denotes the ballistic case). This is a 
clear analogy with the case of Josephson S/F/S junctions with a nonhomogeneous magnetization where 
short- and long-range condensate components are induced in the F-layer. However, in contrast to the charge 
Josephson current in S/F/S junctions, the spin current in AF/F/AF junctions is not constant in space, but 
oscillates in the ballistic F-layer. We also calculate the dependence of ygp on the deviation from the ideal 
nesting in the AF / 1 /AF junctions. The spin current is maximal in the insulating phase of the AF and decreases 
in the metallic phase. It turns to zero at the Neel point when the amplitude of the SDW is zero and changes 
sign for certain values of the detuning parameter. 

PACS numbers: 74.45. +c, 74.50.+r, 75.70.Cn, 74.20.Rp 



I. INTRODUCTION 

Similarity between the magnetic and superconducting 
ordering has been noted long time ago (see Ref. 1 and the 
review articles 2 and 3). In the simplest case of the BCS su- 
perconductivity the order parameter is described by a com- 
plex quantity A = | A| exp(/x) with the amplitude | A| that de- 
termines the Cooper pair density and the phase % related to 
a voltage V via the gauge Josephson equation 



tization M (in the antiferromagnets the vector M is the mag- 
netization in one of the sublattices) . The absolute value of M 
corresponds to the amplitude |A| in superconductors, and 
the angle 6 between the z-axis and the vector M is analo- 
gous, to some extent, to the phase %■ The temporal variation 
of 6 is described by the equation similar^'^ to Eq. (1), 



dtB = -g^B^B- (Mr X Ml)/MrMl 



(3) 



hdtx^qv, 



(1) 



with q - -2\e\, e being the elementary charge. This relation 
shows that, in equilibrium, the gauge-invariant quantity 
equals zero: n - hdtX - 2eV. 

The electrical current j,^ of the condensate is expressed in 
terms of another gauge -invariant quantity — the condensate 
velocity Vs - [hVx - qAlc)l2m as\q - qtisVs- 

One of the most remarkable phenomena arising in su- 
perconducting systems is the Josephson effect. Josephson^^ 
showed that a condensate current 



where g is the gyromagnetic ratio, /iB is the Bohr mag- 
neton, 5B = Br - Bl and Br l is the magnetic induction in 
the respectively right or left F- or AF-layer composing the 
Josephson-like structure. In Eq. (3), the z-axis is chosen to 
be directed along the vector 5B.' 

The Josephson-like effect in F/J/F and AF/J/AF junctions 
was studied in Refs. 4-6. It was found that in these junctions, 
the spin current has the form 



isp = icSp sin(a) , 



(4) 



(2) 



can flow in an S /I / S junction being a periodic function of the 
gauge-invariant phase difference between superconductors 
(p - -(2.mlh) Vgdl, where the limits of the integration R 
and L mean respectively the right and left superconductors 
S separated by an insulating layer I, and jc is the Josephson 
critical current. 

In the case of a magnetic ordering — ferromagnetic (F) or 
antiferromagnetic (AF) — the order parameter is the magne- 



where a = (9r - 0l and 7c,Sp is a constant related to the bar- 
rier transmittance \T\^. 

In the full analogy to the Josephson charge current jj, the 
spin current is nondissipative. The absence of the dissi- 
pation for the spin current and the possibility to use this 
dissipationless current in spintronics has been discussed 
in many papers*^"^^ (see also reviews 12 and 13). Unfortu- 
nately, using the spin currents in spintronics does not seem 
to be straightforward because the analogy between the spin 
and charge currents is not complete, namely, the spin cur- 
rent is not a conserved quantity. ^"^ We will demonstrate the 
validity of this statement considering the proximity effect in 



2 



AF/F structures. In particular, we will calculate the spin cur- 
rent through a Josephson-like AF/F/AF junction. 

It is well known that the Josephson effect in S/F/S junc- 
tions reveals new remarkable properties^^"^^ that are absent 
in the case of more conventional S/I/S or S/N/S junctions. 

First, the condensate functions / penetrating into a dif- 
fusive F-layer decay on a very short range of the order 
^F,d - vDhl^, where D is the classical diffusion coeffi- 
cient of the ferromagnet and is the exchange field; the 
subindex d denotes the diffusive case. This expression 
for is valid if the condition « 1 is fulfilled, where 
T is the momentum relaxation time. Usually, the exchange 
field J6' is much larger than the critical temperature of the 
superconducting transition Tc. In addition, the Cooper 
pairs wave function oscillates in space. It consists of the sin- 
glet component /a, and the triplet component /o with zero 
projection of the total spin on the direction of the magneti- 
zation M. Such oscillations lead to a change of the sign of 
the losephson critical current jc, and the so called ;/r- state is 
realized in the junction. 

Another interesting effect predicted in Ref 20 (see 
also the reviews 18 and 19) and observed in several 
experiments^^"^^ in S/F structures, is related to a new type 
of superconducting correlations. This type of correlations 
called "odd triplet superconductivity" can arise due to the 
proximity effect in the F-layer with a nonhomogeneous 
magnetization. In the case of a uniform magnetization 
in the F-layer, the pair wave function there has both the 
singlet component /sng ~ (i/zj (f)i/'j (0) - {Q)i^\{t)) and the 
triplet component /o ~ (if/\{t)y/i (0) -i- (0)i//t (?)> with zero 
projection of the total spin of the Cooper pair on the di- 
rection of the magnetization. Both components oscillate in 
space and decay on a short length of the order of ^pd- In 
the ballistic case, these components oscillate with the pe- 
riod ~ fF,b = vhl2Je. 

The situation changes drastically if the magnetiza- 
tion in F is not uniform in the nearest vicinity of the 
S/F interface. ^^"^^ In this case, not only the compo- 
nents /sng = /s and /o arise in the F-layer, but also the triplet 
component f\{t) ~ (i/aj (0) -i- i/'j (0)i//f (f)> with nonzero 
projection of the total spin on M is generated. This compo- 
nent is insensitive to the exchange field and penetrates the 
F-layer over a long (compared to ^^f) distance. The penetra- 
tion length does not depend on the exchange field J6' and 
is much larger than ^f- In absence of the spin-orbit interac- 
tion and scattering by magnetic impurities this length may 
become comparable with the penetration length of the su- 
perconducting condensate into the diffusive normal (non- 
magnetic) metal ^n^j - [DhllnT]^'^ . 

The component /i can be called the long-range triplet 
component (LRTC) . It is symmetric in the momentum space 
and is therefore also insensitive to the scattering by ordi- 
nary (nonmagnetic) impurities. In this respect the LRTC 
differs drastically from the usual triplet component /tr 
that describes superconductivity in Sr2Ru04.^^ The latter 
is an antisymmetric function of the momentum p and is 
therefore destroyed by ordinary impurities. The correla- 
tor /i changes sign by permutation of the operators i/'j (f) 



and I/'J (0), as it should be according to the Pauli principle, 
which follows from the fact that the Fermi operators (t) 
and i/Af (0) anticommute at equal times. This means that 
/i (0) ~ (0)i//t (0) + y/^ (0)1/^1 (0)) = 0' i-e- the function /i (f) 
is an odd function of time t or, in other words, it is an odd 
function of the frequency w in the Fourier representation. 
That is why this component is called "odd triplet compo- 
nent". 

Recent experiments unambiguously confirmed the exis- 
tence of the LRTC in the S/F structures with a nonhomoge- 
neous magnetization.^^"^^ 

In this paper we show that the analogy between 
the superconducting and magnetic order parameters in 
Josephson-like junctions is deeper than it was indicated 
previously. It turns out that in the SDW/F/SDW junction 
with a ballistic F-layer the "short"- and long-range AF or- 
der parameters are induced in the F-layer in an AF/F struc- 
ture (we use the quotation marks for "short" -range, to in- 
dicate that in the considered ballistic case this component 
does not decay on a short range, but rather oscillates with 
a short period). The component of the AF order parameter 
with the M vector (in each sublattice) being parallel to the 
magnetization vector in the F-layer Mf, which is supposed 
to be oriented along the z-axis, penetrates into the F-layer 
and oscillates with a short period of the order of ^pb- On 
the other hand, the SDW with the M vector perpendicular 
to Mf penetrates the F-layer over a long distance of the or- 
der of ^N_b and decays monotonously. This part of the SDW 
order parameter may be called the long-range component 
induced in the ferromagnetic layer. 

At the same time, in spite of the long-range penetration 
of the AF correlations into the ferromagnet F, the spin cur- 
rent that arises at the AF/F interface is not constant in space. 
It oscillates in the F region with a short period. The study 
of the spin current in various heterostructures is important 
not only from the point of view of fundamental physics but 
also from the point of view of possible applications of these 
structures in spintronics.'^^ 

We do not analyze the case of the Bose-Einstein con- 
densation in the magnon system at a high magnon density 
which can be realized via the external pumping. In this case, 
a Josephson-like effect may also arise, but it has different 
properties (see the review 14 and references therein). 



II. MODEL 

In order to describe the AF order we adopt the model 
developed for superconducting Fe-based pnictides.^^"^^ In 
this model the band structure is assumed to consist of two 
bands (electron and hole) with perfect or almost perfect 
nesting, which is taken into account by a certain param- 
eter Sjd. Due to this, a spin density wave (SDW) corre- 
sponding to an antiferromagnetic order parameter arises in 
these materials. On the basis of this model, which describes 
well many properties of this new class of high- Tc supercon- 
ductors (see review articles 36-41 and references therein), 
the generalized Eilenberger equations for the quasiclassi- 
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FIG. 1. (Color online.) The considered setup. The two antiferro- 
magnets SDW are separated by a thin ferromagnetic layer of thick- 
ness 2L. The orientation of the magnetization vectors of the SDW 
in the leads lying in the y-z-plane is shown in red. 

cal Green's functions g have been derived in the previous 
paper'*^ by the authors (see also a recent paper^^ where the 
definition of the Green's functions differs from ours). The 
quasiclassical formalism is especially suitable for tackling 
nonhomogeneous problems and the matrix functions g al- 
low one to find all the necessary observable quantities, such 
as the charge and spin currents, the density- of- states, etc. 
The technique of quasiclassical Green's functions is pre- 
sented, e.g., in Refs. 43-45. 

III. SDW/F/SDW "JOSEPHSON" JUNCTION 

A. "Short"- and Long-range Components 

First, we consider a Josephson-like AF/F/AF junction 
composed of two leads and a thin ferromagnetic layer be- 
tween the leads. In the leads a spin density wave exists, 
which is a particular example of the AF order and the layer 
between these materials is a ferromagnet in the ballistic 
regime (see Fig. 1). For simplicity, we assume that both the 
SDW and the F parts of the structure consist of similar two- 
band materials and differ from each other by different in- 
teraction constants and doping levels. The AF/F interface 
penetrability is supposed to be small and thus one can ne- 
glect the effect of the F-layer on the SDW leads. The Green's 
functions in the right (left) leads are given by the formula 
that can be obtained from Eq. (6.4) of Ref 42 by setting the 
superconducting order parameter to zero, A = 0) 

1 , , 

gR(L) = -— {<y„T3 + WM0[pi-T2-O-3COS(a/2) (5) 
&M 

+ P2 -Ti -(72 sin(a/2)]}. 

Here, (y„ - nT{2n + 1) is the Matsubara frequency, Wmo is 
the amplitude of the SDW, <Sm - ^J(t^i + W^ and p, • fj ■ (T^ 
is the tensor product of Pauli matrices in the band, Gor'kov- 
Nambu and spin space, respectively. The magnetization of 
the SDW lies in the y-z-plane at the angle +a/2 relative to 
the z-axis, with the projection on the y-axis +Wjvfosin(a;/2) 
in the left (right) leads. For simplicity, we suppose the ideal 



nesting and set dj^-O. In the F-layer, the quasiclassical 
Green's functions obey the equation 

n^vd^g+[A,g]^Q, (6) 

where nx- Pxlp and v is the modulus of the Fermi veloc- 
ity. The matrix A is defined as A = f 3 • {(UnPa ■ + i-^Po • ^'3) 
with the exchange field in the ferromagnet Jf?. We assume 
that there is no impurity scattering in the F-layer (ballistic 

case). The impurity scattering leads to a suppression of the 
SDWamplitude.51'52 

Equation (6) is supplemented by the boundary condition 
which we write in the form 

(g(«x) - gi-rix)) = ±sgn(«;,) [g,#' ■ gR(L) ■ . (7) 

The transmission coefficient 3' is & matrix, which is as- 
sumed to be small. This boundary condition generalizes 
the widely used boundary condition derived by Zaitsev*^ 
to the case of a two-band metal with an SDW. We con- 
sider a simplified model for the AF/F interface and assume 
that there are no interband transitions. In this case, the 
matrix 5" has the form 3' - (^p3 -1- ^po) -f -CTq, where 
5o,3 - {Ti + 72)12 and T12 are the matrix tunneling elements 
for transitions from the band 1 resp. 2 of the F-layer to the 
same bands 1 resp. 2 of the SDW leads (see the Appendix A). 
If these elements are equal, i.e., T\-T2= Ttun, then 3o-0 
and ^3 = Ttun. 

The boundary conditions that describe real materials are 
more complicated. They must include the transitions be- 
tween different bands and this process may be accompa- 
nied by spin flips. However, we are not interested in exact 
calculations of the amplitude of the SDW penetrating into 
the ferromagnetic layer, but rather in a qualitative effect — 
the appearance of "short"- and long-range components. 
The existence of these components does not depend on the 
exact form of the boundary conditions. 

In the lowest order of approximation in the transmission 
coefficient 3' (no proximity effect) , the matrix g in the ferro- 
magnet has the form g = gop - sgn((y)po ■ f 3 • ctq. Therefore, 
the boundary condition can be written in the form 

iginx) - gi-rix)) = ±sgn(n;c)^eff[*3'/R(L)] , (8) 

where the matrix /r(l) equals 

/r(L) = — — {pi-f2-(T3Cos(a/2) + p2-fi-(t2sin(a/2)} (9) 
6m 

and 

^eff = - (-^O^ - -^3^) ^TiT2 (10) 

is an effective transmission coefficient for the SDW. The 
boundary condition (8) defines the spin current through 
the interface. One can see that if one of the tunneling ele- 
ments {Ti or T2) is zero, then the spin current turns to zero. 

We need to find a solution of Eq. (6) supplemented with 
the boundary condition (8). In the considered case of 
a low penetrability of the SDW/F interface, the matrix g 
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in the F-film can be represented in the form g - goF + Sg, 
where dg is a matrix with a small "amplitude" which de- 
scribes the SDW induced in the F-film due to the proxim- 
ity effect. As usual, we represent g as a sum of a sym- 
metric in momentum part and the antisymmetric one, i.e., 
Sg - s + rixa. We substitute this expression into Eq. (6) and 
split it into two equations for the symmetric and antisym- 
metric parts 



vdxS+ [a, a] - 0, 
n\vdxa+ [A, 5] =0. 



(11) 
(12) 



Excluding the matrix a, we obtain one equation for the ma- 
trix s 



-v^nldls+[A,[A,s]\^Q. 



(13) 



For the matrix a we obtain a similar equation. The commu- 
tator [ A , [ A , 1] ] is easily calculated resulting in 

[A,[A,S]] =2{(W^-^^)S-H2i^W((T3-S-HS-(T3) (14) 

Here, we used the property 5~fi2- It is seen from the 
boundary conditions (8-9), that this matrix anticommutes 
with 73. It is convenient to represent the matrix 5 as a 
sum of two matrices 5 = 5|| -1- 5^, where s\\ = sq • <to -1- S3 • 173 
and = 52 ■ (72. We will see that the matrix s\\ describes a 
"short"-range component of the SDW penetrating into the 
F-layer, while sx_ determines the long-range component, 
which is perpendicular to the orientation of the magnetiza- 
tion in the ferromagnet. From Eq. (13) we obtain the equa- 
tions for 5|| and sj_ 

-v^nld\s\\+A{(x) + \jea2fs\\^Q, (15) 
-v^nld\s^ + A(jy^s^^Q. (16) 

In the first order in the transmission coefficient ^T^^ the 
solutions of these equations obeying the boundary condi- 
tions are 

S|| - 2mSa,3~^ffPi ■ ■ [(73 ■S?(x) -I- i&o ■J'{x)\ ■ cos(q;/2) , 



S]_ - 2mSa,3'^ffP2 • f 1 • (T2 • (x) • sin(a/2) , 



(17) 
(18) 



where Sa, = sgn(<i)) and m - WmqI^m is the dimensionless 
amplitude of the SDW induced in the F-layer and 



Sg(x) = !ft{^(x)}, 



J?(x) = !3{^(x)} 



and 



cosh(0+x/L) 
sinh(0+) 



cosh(0ft,) 



(19) 



(20) 



with = KojL and 9+ - k+L. The here introduced wave vec- 
tors )c+ = 2 (|(y| -i-i^sgn(w))/L'|nj;| and - 2\(i)\l v\nx\ de- 
termine the characteristic lengths which describe the pen- 
etration of the "short"- and long-range components of the 



SDW into the F-region, respectively. The same wave vectors 
characterize the penetration of the short-range component, 
i.e., the singlet and triplet component with zero projection 
of the total spin onto the z-axis, and the long-range triplet 
component of the superconducting wave function into the 
F-region of an S/F bUayer in the ballistic case.^''"^^ 

The component of the SDW Eq. (15), which can be 
called the "short" -range component, oscillates at a given 
angle x - arccos(nj;) rather fast in space with a period of the 
order of ^p.b (it is assumed that ^» T). Being averaged 
over the angle, it decays from the interface in a power-law 
fashion. The long-range component s_i, Eq. (16), penetrates 
into the ferromagnet over a long length ~ ^^ b, which does 
not depend on the exchange field J^. Thus, there is an anal- 
ogy with an S/F/S Josephson junction where the ferromag- 
net has an inhomogeneous magnetization.^^'^^ In the bal- 
listic case, the "short" -range component of the condensate 
wave function oscillates in space with a period ~ fp.b/'^"^^ 
whereas the long-range component penetrates the ferro- 
magnet over the length ~ ^^N,b- We use quotation marks for 
the "short" -range component because, in the considered 
ballistic case, this component does not decay exponentially 
in the ferromagnetic layer, but oscillates and decays in a 
power-law fashion after averaging over the angles. In the 
diffusive limit in the S/F/S junctions, this component de- 
cays exponentially on a short distance of the order 

Note that the "short"-range component S|| consists of the 
odd and even in w parts (correspondingly, the first and the 
second terms in Eq. (17)), whereas the long-range compo- 
nent s±_ is an odd function of the frequency w. The same 
dependence has an anomalous (Gor'kov's) Green's function 
in an S/F/S system with a nonhomogeneous magnetiza- 
tion M.^^'^^ In the last case, the short-range part consists 
of the singlet component and triplet component with zero 
projection of the total spin onto the M vector, and the long- 
range part consists of the triplet component with a non-zero 
projection of the total spin onto the M vector. 

In Fig. 2 a) we plot the coordinate dependence of 
the "short"- (the functions S^{x) and J^(x)) and the 
long-range (the function Sg'(x)) components. The former 
stems from the part of the matrix g corresponding to the 
magnetization component of the SDW M parallel to the 
magnetization in the ferromagnet Mf||z, i.e., the second 
term on the right in Eq. (5). The long-range compo- 
nent stems from the SDW component perpendicular to the 
Mp vector and turns to zero at a = 0. 

The components of a that determine the spin cur- 
rent ]'^p^\x] can be presented in the same form as the 
s\i± matrix, i.e., a - a\i + d±. The "longitudinal" — dp and 
the "transverse" — d±_ , parts are related with S|| ^ through the 
simple formulas (see Eq. (12)) 



vnldxdii = 2S(uf3 • (|w|(to -1- i^&i) ■ s\\ , (21) 
vnldxd± = 2soj\cD\f3-ao-s_L. (22) 
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FIG. 2. (Color online.) The coordinate dependence of the "short"- 
range resp. the long-range components of the Green's function in 
the F-layer. 

a) the symmetric in momentum part of g — 5?(x) and J'(x) 
resp. "^[x) 

b) the antisymmetric in momentum part of g — ^ix) and J^(x) 
resp. "^[x) 



From Eqs. (15-22) we find the matrix a: 
■fa • = -2m— ^pi - fa ■ [cts -^(x) + i(7o • J^(x)] ■cos(a/2), 



\nx\ 



a 



eff 

T3 ■ fli = -2m p2 ■ Ti ■a2-^{x) ■ sin(Q:/2) , 

l"xl 



where 



J^(X) = !3{/(X)}, 



(23) 
(24) 

(25) 



and 



/{x) = 



sinh(0+x/L) 
sinh(0+) 



cosh(0,,,x/L) 

<^(x) = . (26) 

cosh(0(j) 



These functions are plotted in Fig. 2 b) and contribute to 
the "short"- and the long-range components of the Green's 
function in the F-layer, more exactly — to the odd in rix part 
of it. 

As it follows from the expression for the spin current (see 
Eq. (3.13) in Ref. 42)), the obtained formulas for ai±_ lead to 
zero spin current inside the F-film. This result is quite natu- 
ral because the spin Josephson-like current in SDW/F/SDW 
junction should be proportional to the penetration proba- 



bilities of both interfaces, i.e., ]'^p^\x] 

symmetric case under consideration 3'^^ 
In order to obtain a finite spin current, we have to calcu- 
late the matrix a in the next order in the transmission co- 
efficient 3'^2- '-^^ easily find the corrections da\\^±_ from 
the normalization condition 



. ar2 _ ar-l 
eff (L) ~ ^ eff(R) ~ ^ eff ' 



= 1. 



(27) 



where at - a + 8a and St-s + ds. 

Therefore, for 6a we obtain the equation 

2goF -80. + a-s+s- a-Q, 



(29) 



where goF - sgn(w)po ■ f 3 ■ (Tq. We are interested in the com- 
ponent 5fl3 of the matrix Sa, which is proportional to the 
matrbc 73, because only this component contributes to the 
spin current. This component commutes with the ma- 
trix goF- 

Substituting expressions (17-26) into Eq. (29), we obtain 

for 5(23 



('"■^eff)^ „ I cosh [(0+ - B^)xlL] 



Sa-i- 2is„ 



\nx\ ( sinh(0+)cosh(0„) 
One can write Eq. (30) in the form 

cos(0^) 



sin(a)p3-f3-(Ti. 

(30) 



da'i — 2\Si 



1 tanh(t^(y) „ 

\nx\ cosh2(0^)-cos2(0^) 

X cos{6 jexl L) sin(a;)j03 • 73 -cti , 



(31) 



where Qje- J€LI v\nx\. 

Now, the spin current can be easily calculated. 



B. Spin Current 

We see that the SDW/F/SDW junctions are very simi- 
lar to the S/F/S Josephson junctions. In analogy to the 
"short"- and long-range components of the superconduct- 
ing condensate penetrating the ferromagnet, there are sim- 
ilar "short"- and long-range components of the SDW pene- 
trating the ferromagnet. 

Nevertheless, there is an essential difference between the 
two systems. The Josephson charge current in an S/F/S 
junction is constant in space in the ferromagnet. It is de- 
termined by the long-range component, provided the mag- 
netization M(x) in the F-layer is not uniform and the thick- 
ness of this layer is not too small (dp » ^f)- 

In contrast, the spin current in an SDW/F/SDW junction 
is not constant in the F-layer but oscillates with a period of 
the order of (see Fig. 3). It is constant in space only in the 
absence of the exchange field, i.e., in SDW/ N/ SDW junction. 

The spin current through the SDW/F/SDW junction can 
easily be calculated using the formula for jsp (see Eq. (3.13) 
of Ref. 42) 



X2nT)vv 



The notation (x) means the spin current density with 
the spin projection x (resp. y] flowing perpendicular to the 
S/F interface and the angle brackets denote the angle aver- 
aging. 

Substituting Eq. (30) into Eq. (32) we obtain 



|"Bl^Tr<(p3 -73 -(71,2 •«)«?>. (32) 



As it follows from this equation 



at ■ (goF + St) + (goF + St)-dt^O, 



(28) 



ispW = ic':spW-sin(a:), 



(33) 
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FIG. 3. (Color online.) The coordinate dependence of the spin cur- 
rent in the F-layer for different values of the strength of the mag- 
netic field h = J£; T,„ is the SDW transition temperature. 

a) (5-eff)2(X 5(l-|«xl), 

b) (^eff)^ =TB = const. 



where 



jc,SpM ^ fiBVV[2TlT) 



(34) 



xE 



I nx\m'^ 3' *^^tanh{0 a,) cos(0^) cos(d jexIL) 



cosh^(0w)-cos2(0^) 



One can see that the spin current oscillates in space in the 
F-layer. It is also seen that in the case of a normal (nonmag- 
netic) metal, i.e., at = 0, the spin current is constant in 
space. 

As follows from Eq. (34), the spin current at the interfaces 
it is equal to 



|wUm2^Jjtanh(0a))cos2(0^) \ 
cosh2(0^)-cos2(0^) /■ 



(35) 



Using Eq. (34) we plot the coordinate dependence of the 
critical spin current with the projection of the spin onto the 
X-axis in Fig. 3 for different values of The critical spin 
current oscillates with a period of the order of ^pb- On the 
left panel we choose (.5eff] oc S{\-\nx\) so that there is no 
dependence of the oscillation amplitude on x; on the right 
panel we choose (^ff)^ -Jb- const., so that the amplitude 
depends on x in a power-law fashion, which is a conse- 
quence of the integration over the angles. 

Thus, we see that the SDW order parameter penetrat- 
ing into the ferromagnet contains two components — a 
long-range component with M vector perpendicular to the 
magnetization in the ferromagnet Mp, and a "short"-range 
component with M vector parallel to the z-axis. A simi- 
lar effect arises in an S/F/S junction with the nonhomoge- 
neous magnetization, where the long-range component is 
the odd triplet (S ^ 0) wave function of the Cooper pairs, 
and the short-range component corresponds to the sin- 
glet and S-Q triplet wave function. ^^'^^ However, the spin 
current jsp is not analogous to the charge Josephson cur- 
rent j] — it is not constant in space, but oscillates with a pe- 
riod ~ ^pb- 



IV. SDW/I/SDW "JOSEPHSON" JUNCTION 
A. Ferrell-Prange Equation 

Now, we consider a junction formed by two materials 
with the SDW. These can be iron-based pnictides separated 
by an insulating, normal metal or ferromagnetic layer. In all 
these cases the current through the junction is given by 



isp = ic':sp-sin(«). 



(36) 



where the critical spin current depends on the type of the 
barrier. For example, for SDW/N/SDW junction it is equal 
to 



ic"sp = 2/XBJ^v(2;rr)M^^oX 



<f^sinh(20^) 



(37) 



Here, the superscript x denotes the component of the spin 
current vector having the spin projection on the x-axis only. 
For the case of SDW/I/SDW junction we easily obtain 



/c"sp = 2MBi^v(27rr)<o^ 



(38) 



Introducing f = {\n\x3'^^ we can write this equation in the 
case of ideal nesting as 



Jc.Sp 



- - tanh ■ 



M 



(39) 



27tT 



which corresponds to the Ambegaokar-Baratoff formula for 
the critical Josephson current with A replaced by Wmo and 
the "interface resistance per unit area" R is related to the 
effective transmission coefficient — R oc 1/f . 

Equation (36) for the spin "Josephson" current in an 
AF/J/AF junction has been obtained in Ref. 6 by another 
method. 

Now, we derive an equation similar to the FerreU- 
Prange*^ equation. We consider a junction with SDW leads 
in the form of thin films. The thickness d of the films must 
be less than a characteristic length, over which the SDW vec- 
tor M restores its preferable direction. This direction is de- 
termined by anisotropy effects. 

The idea of the derivation is the same as in Ref. 48. Con- 
sidering the region near the boundary (cf. Fig. 4) we use the 
continuity equation for the spin M''^^ 



afM'^) = -divjsp, 



(40) 



which can be derived easily from Eq. (3.5) of Ref. 42. Using 
the Green's functions g unperturbed by the proximity effect 
(see Eq. (5)), one can show that M'^' - 0. From the continu- 
ity equation we get 



(41) 



The formula for the spin current j^^ in an SDW lead with 
the magnetization orientation varying in space can be easily 
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where the size of the "fluxon" is given by Iq - y/ die. The 
mutual orientation changes along the z-axis from to 2n 
and has the value n in the core of the "fluxon", cf Fig. 5. 
In the case of the Josephson fluxon one has a circulating 
charge current — in our case, the spin current circulates in 
the "fluxon". 



FIG. 4. (Color online.) A scheme for derivation of the Ferrell- 
Prange equation. 



^ ._t J ..J^.. 



— CO 



^ 7 1 "^ "V 



B. Finite Detuning 

Up to now, we assumed the ideal nesting, 5/i - 0. In the 
case Sj^^O we find the quasiclassical Green's function in the 
SDW system by inverting the matrix G in Eq. (3.3) of Ref. 42. 
It turns out that g is a linear combination of six "basis ma- 
trices" 

gL/R - g030 ± gl22 + gl23 ± g212 + g213 + gSOO . (47) 



FIG. 5. (Color online.) Dependence of the magnetization direction 
in the "fluxon" with z. The mutual orientation in the leads changes 
from to 2ji. 



where the subscripts L and R denote the left and right SDW 
systems, respectively, corresponding to the setup depicted 
in Fig. I with the F-layer replaced by a very thin insulating 
barrier. The six functions gijk are 



obtained from Eq. (6) in the same way as the Meissner cur- 
rent is found in a clean superconductor (see, e.g., Ref. 45). It 
has the form (see Appendix B) 



7sp(^)|. 



1 da 

6 az 



M 



(42) 



This equation formally coincides with the corresponding 
formula for the supercurrent in a clean superconductor if 
one replaces Wmo ^ A, % ^^d /i^ 2e, where % is the 
phase in the superconductor. 

Substituting Eq. (36) into Eq. (41) we arrive at the final re- 
sult 



dlcdia - sin(a) , 



(43) 



where Ic is a characteristic length inversely proportional to 
the critical current jc.sp- It equals 



eff 



<?^sinh(20^) 



0)) I I 



(44) 



g030 = ll+ ■ SR + } -K5p ■ 3 {x+]] • po • f 3 ■ (To , (48) 

gi22 = i| 1+ T^IVmo • 9 {x+} sin(a/2) ■ pi ■ f 2 • (T2 , (49) 

gi23 = I1+ IVmo ■ 5ft {1+ } cos(a/2) ■ pi ■ f 2 ■ (73 , (50) 

g2i2 = ll+r^WMo-5ft{l+}sin(a/2)-p2-fi-(T2, (51) 

g2i3 = \\x+\~^Wmo ■ 9 {1+} cos(a/2) ■ p2 • f 1 ■ (T3 , (52) 

g300 = i|l+ (w • 9 {x+} - • SR • p3 ■ f ■ (To , (53) 



where X+^\l W^mo + ^'^ + i^^^'^- 

The detuning parameter Sfi can be written in the form 
S fi - fiQ + fi(f, cos(,2(p) with fiQ describing the relative size 
mismatch of the "hole" and "electron" pockets and /i^ de- 
scribes the ellipticity of the "electron" pocket, while the 
"hole" pocket is assumed to be a circle.^^'^^'^* 

From the quasiclassical Green's function, Eq. (47), we can 
calculate the dependence of various quantities of interest 
on Sfi. 

We begin with the density of states (DoS). The DoS is 
given by the formula 



in the case of SDW/N/SDW junctions and 



hv 3 



(45) 



for SDW/I/SDW junctions. 

Equation (43) is the analog of the Ferrell-Prange equation 
for a junction composed of two materials with the SDW and 
separated by some barrier (a thin insulating or normal metal 
layer). 

One can easily find a solution similar to the one describ- 
ing a fluxon in the losephson tunnel junctions. It has the 
form 



a - 4arctan [exp(z/ Iq)] 



(46) 



v(e) = i$R(Tr(po-f3-(To-g(£-))), 



(54) 



where the angle brackets mean the averaging over the mo- 
mentum directions 



(55) 



Inserting the expression for the quasiclassical Green's 
function into Eq. (54) we obtain 



I 

2 



- 5fi 



(C-K5p)2 



MO 



(56) 



8 





FIG. 6. 



(Color online." 
■■ 1.07. 



DoS for different values of s = fi^ljio at 




FIG. 7. (Color online.) WMoiH(p) at /io/Js = 1.07; black dashed 
curve shows the dependence of DoS on /j^ at zero energy. 



This result was obtained earlier in Refs. 33 and 34. 

We plot the DoS for different values of the ratio s = /i^/po 
in Fig. 6 choosing a constant /io - LOTTs, where Ts is the 
SDW order transition temperature. Note that there is a 
range of the values of 5 where an energy gap appears in 
the DoS, indicating an insulating phase of the system. We 
can trace the transition to the metal considering the value 
of the DoS at zero energy, v(0), as a function of the param- 
eter fifj,. In Fig. 7 we plot the dependence of the order pa- 
rameter Wmo on and, in the same plot, show the depen- 
dence v(0) on /i0 (dashed line). One can observe a rather 
steep increase of v(0) in the region around 0.7rs. 

The dependence of the SDW order parameter Wmo on /x^ 
is calculated from the self-consistency equation, Eq. (3.19) 
of Ref. 42, which in addition has to be averaged over the mo- 
mentum directions. It reduces to the form (see also Ref. 33) 



ln(r/rs) = 27rr^5R(x- 

w>0 



- 0) 



(57) 



It is of special interest^ to calculate the spin current 
through the SDW/ 1 /SDW system. Using the expression 
Eq. (47) for the quasiclassical Green's functions in the left 



FIG. 8. (Color online.) Spin current on /io resp. /i^. 



and right leads, and inserting it into the expression for the 
spin current, Eq. (3.13) of Ref. 42,we obtain the familiar form 
of the Josephson-like current 



JSp = 7c,Sp-sin(a;) 



(58) 



for the spin current through the SDW/ SDW interface which 
has only the spin projection on the x-axis. The critical coef- 
ficient 7c,sp is dependent on 5/i 



7c,Sp = 2^bVVWIjq[271T)( iHxI^eff E 
\ aj>0 

(59) 

In Fig. 8 we show the dependence of the critical spin cur- 
rent on the two parameters fiQ and fi^f,, choosing in each case 
the other parameter constant — when plotted as a function 
of /io, we set /i^ - 1.26rs; plotting the spin current as a func- 
tion of /i0 we set /io - 1.07 Tg. 

One can see from Fig. 8 and Fig. 7 that in the insulating 
phase the spin current is not zero. In contrary, it is at its 
maximum. This addresses the question about the possibil- 
ity of the appearance of the spin current in an insulator (see 
discussion on p. 12 of Ref. 3). 

Analyzing the expression (59) for the critical spin current 
we see that there is a possibility of a sign change. At small 
temperatures, the spin current changes sign if, e.g., the con- 
dition Wmo /"o - M0 holds with /x^ < /io at the same time. 

We note the similarity of the self-consistency equa- 
tion (57) with the self- consistency equations presented in 
Refs. 49 and 50. Setting /i^ = we can map the problem 
of finding the dependence of Wmo on /io onto the prob- 
lem of finding the dependence of the superconducting or- 
der parameter on a strong exchange field resulting in the 
so-called Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) state, 
which is characterized by the spatial dependence of A. Note 
that the possibility of the LOFF-like state in pnictides was 
noted earlier by Gor'kov and Teitel'baum.^^ 



V. DISCUSSION 

We have considered Josephson-like junctions of the 
SDW/F/SDW and SDW/I/SDW types and calculated the 
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spin current jsp in these systems. In both the cases the de- 
pendence of 7sp on the angle a between the mutual orien- 
tations of the magnetization vectors M in the right and left 
leads is given by Eq. (36) . The critical current density jsp.c is 
determined by the interface transparencies. 

We emphasize the analogy between the SDW/F/SDW 
junctions and the Josephson S/F/S junctions with a non- 
homogeneous magnetization. In both the systems there 
are "short"- and long-range components of the order pa- 
rameters that penetrate the ferromagnet. In the ballistic 
systems under consideration the "short" -range SDW com- 
ponent penetrating the F-layer oscillates with a period of 
the order of ^p,b - hvlJf. It stems from the SDW compo- 
nent parallel to the magnetization vector Mp in the ferro- 
magnet. The SDW component normal to the Mp vector 
penetrates the ferromagnet on a long distance of the or- 
der of (J^N.b - hvllnT. The penetration depth of this com- 
ponent does not depend on the exchange field Jf?, and, in 
this sense, this component is analogous to the odd triplet 
long-range component in S/F/S junctions with a nonhomo- 
geneous magnetization. 

However, there is an essential difference between the spin 
current in SDW/F/SDW junctions 7sp(^) and the Joseph- 
son charge current j] in S/F/S junctions. Whereas the 
Josephson current j], as it should be, does not depend on x 
(divjj - dxj] - 0), the spin current 7sp(^) oscillates in space 
with a period of the order ^p.b- This behavior is a conse- 
quence of the fact that the spin current arises due to the in- 
terference between the "short"- and long-range component 
of the SDW in the ferromagnet. 

The spin current in SDW/I/SDW junctions depends on 
the angle a in the same way as in the SDW/F/SDW junc- 
tions, i.e., 7sp ~ sin(a). We have calculated the dependence 
of 7sp in an SDW/ 1 /SDW junction on the parameter 5/i char- 
acterizing the deviation from the ideal nesting. At 6 p. - 0, 
an energy gap Wmo opens in the excitation spectrum in the 
leads with the SDW (see Fig. 7). In a certain region of the pa- 
rameter the energy gap disappears and the leads of the 
junction become metals (or semimetals). The spin current 
exists in both the cases; it turns to zero at the Neel tempera- 
ture when Wmo 0. 

Furthermore, we have derived the Ferrell-Prange equa- 
tion for the Josephson-like junctions with the SDW and 
found a solution describing a localized spin current distri- 
bution analogously to the fluxon in a tunnel S/I/S Joseph- 
son junction. The characteristic length of the "fluxon", 
lo = \/dlc, is determined by the barrier transmittance and 
the thickness d of the leads with the SDW. In the case of bulk 
leads the thickness d should be replaced by a characteris- 
tic length over which the vector M of the SDW restores its 
favorable direction in the bulk, which is determined by the 
anisotropy effects. 
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Appendix A: Model Boundary Condition 

Here, we generalize the boundary conditions for quasi- 
classical Green's functions to the case of a two-band ma- 
terial with the SDW. It is not in the scope of the paper to 
present a rigorous derivation of the boundary conditions 
with account for all the possible processes accompanying 
the passage of electrons through the interface. We ignore 
the interband transitions and also the spin-flip processes. 
In addition, we restrict ourselves to the case of a low trans- 
mittance coefficient. Thus, the derived boundary condi- 
tions, to some extent, are phenomenological. However, in 
a single-band case they coincide with the boundary condi- 
tions derived by Zaitsev^^ and in case of two-band material 
they are compatible with the boundary conditions used in 
Ref. 53. Note also that the exact form of the boundary con- 
ditions is not important for our qualitative conclusions such 
as the appearance of the "short" - and long-range SDW com- 
ponents in the ferromagnet. 

In the case of weak penetration and single-band materi- 
als, the boundary conditions for the Green's function g have 
the form of Eq. (7), in which the matrix 3' should be re- 
placed by a scalar 5" (see Ref. 46). The Green's function on 
the left-hand side of this equation determines the charge 
current through the interface. It equals the commutator 
multiplied by 3'^ and consisting of the symmetric in mo- 
mentum space Green's functions. This part of Eq. (7), which 
of course also determines the charge current, can be ob- 
tained by using the tunneling Hamiltonian method (see, e.g, 
Ref. 54). In terms of the operators Ci(j) and C^^.^ introduced 
in Ref. 42 the tunneling Hamiltonian can be written as fol- 
lows 

E (^C/Cr-HCcj , (Al) 
k,s,s' 

where 3' - [S^ps + 3opo) • fs • ctq and the coefficients i5o_3 
are defined in Eq. (7) in terms of the matrix elements Ti 2 
describing the electron transition through the interface in 
each band. The Green's function G satisfies the equation 
(see, e.g, Ref. 54) 

Gi = Go,i+Go,i-Zt-Gi, (A2) 

where l.j -3~Gi3'. This means that the self- energy part 
has formally the same form as in the case of the impu- 
rity scattering, but the Green's function in l.j describes an- 
other (right) lead in the equation for the Green's function Q 
of the left lead. The charge current through the interface is 
determined by the trace of a commutator 

Tr{p3 ■ f 3 • CTo ■ [It , Gi] } = Tr{p3 • f 3 • ctq ■ [^Gr^" , Gi]} ■ 

(A3) 

The expression in the curly brackets expressed in terms of 
the quasiclassical Green's functions appears in the right- 
hand side of Eq. (7). 



10 



Appendix B: Spin Current in a Bulk Material with the SDW 

In order to obtain the expression for the spin current in 
the left (right) lead flowing parallel to the interface, we as- 
sume that the angle a varies in space slowly so that the 
change of a on a short coherence length ^^sdw - hu/WMo is 
small. The characteristic length of the a variation is in- 
versely proportional to the critical current j^^^ (see below) 
and is much larger than f sow- If the angle a is spatially con- 
stant, then the Green's function g = ga obeys the equation 

[A,g„]=0. (BI) 

The solution of this equation is given by Eq. (5) and may be 
represented in the form 

ga^AISM^SgoS\ (B2) 

where the matrix go = {w^fa -i- (Wmo/<S'm)Pi •■r2 -cts} does 
not depend on a. The unitary matrix S is defined as 
S = cos(a/2) -I- isin(Q;/2)p3 ■ fs -(Ti. If the angle a(z) slowly 



depends on the coordinate along the interface, we have 
to add a correction dga to the matrix ga which depends 
on rix — 5ga - rixOa- The matrix a satisfies the equation 

v{S^-d,S-go + go- d,S^ ■ S) + 2SMgo ■ flo = , (B3) 

where the matrices and uq are related in the same way as 
matrices go and ga, i.e., Ua = SoqS^. Here, we used 

go-ao + ao-go = 0, (B4) 

which follows from the normalization condition Eq. (27). 
One easily finds 

S^dJ^[il2)d,ap3-f3-(Ji- (B5) 

Combining Eqs. (B3), (B5) and the normalization condi- 
tion go • go = 1, we find 

ao-- v(il2)dzaWM0<S^{- iopo ■73-^2 + WmoPs • fa • (ti] , 

(B6) 

and using Eq. (32), we obtain the Eq. (42) for j^^iz). 
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